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Abstract 

We propose a new exactly solvable model of strongly correlated electrons. The 
model is based on a d-p model of the Cu02 plane with infinitely large repulsive inter- 
actions on Cu-sites, and it contains additional correlated- hopping, pair-hopping and 
charge-charge interactions of electrons. For even numbers of electrons less than or equal 
to 2/3-filling, we construct the exact ground states of the model, all of which have the 
same energy and each of which is the unique ground state for a fixed electron number. 
It is shown that these ground states are the resonating-valence-bond states which are 
also regarded as condensed states in which all electrons are in a single two-electron 
state. We also show that the ground states exhibit off-diagonal long-range order. 
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1 Introduction 



Considerable attention has been paid to strongly correlated-electron systems in the study 
of high-T c superconductivity since Anderson pQ proposed that its origin may be attributed 
to magnetic interactions induced by strong Coulomb repulsion. In Anderson's scenario, 
a Mott insulating phase of an undoped cuprate is assumed to be the resonating-valence- 
bond (RVB) state, composed of electron spin-singlet pairs, and these singlet pairs become 
charged superconducting pairs when the Mott insulator is doped sufficiently Although 
many attempts have been made on lattice models of electrons with interactions, such as the 
Hubbard model and the t-J model, there is no definite answer to the question of whether this 
scenario is realized in these concrete microscopic models. These models are usually difficult 
to analyze theoretically, and exact results are limited except for one dimensional case where 
Bethe ansatz is available. 

One of the exact results in two or more dimensions is obtained by Brandt and Giesekus l2 j , 
who constructed the exact ground state of the Hubbard model with infinitely large repulsive 
on-site interactions on a <i-dimensional decorated hypercubic lattice. Tasaki |3] generalized 
the Brandt-Giesekus model to a class of Hubbard models on certain cell structures. He 
pointed out that the exact ground states in this class of models are regarded as RVB states. 
He also obtained the singlet pair correlation function on the tree lattice and discussed a 
possibility of long-range order associated with the singlet pairs and superconductivity. Bares 
and Lee [3] and Yamanaka et al. |3] treated one dimensional versions of the models and 
obtained various correlation functions, including the singlet pair correlation function. They 
found that all the correlation functions decay exponentially with distance and concluded that 
the system is in an insulating phase. There are no conclusive results about the occurrence 
(or absence) of superconductivity in the models in two or more dimensions. 

For extended Hubbard models including nearest-neighbor interactions in addition to the 
on-site interactions, there are some exact results relevant to superconductivity. These results 
are based on the so-called 77-pairing mechanism originally proposed by Yang [H] . By using the 
77-pairing mechanism Yang constructed eigenstates of the Hubbard model on the hypercubic 
lattice and proved that these states, in which all electrons form the 77-pairs, have off-diagonal 
long-range order and thus are superconducting. Unfortunately, these eigenstates were shown 
not to be the ground states. But, Essler et al. [7] showed that the ^-pairing states become 
the exact ground states of an extended Hubbard model with attractive on-site interactions. 
They also showed that the model with moderately repulsive on-site interactions has ground 
states in which part of electrons forms the 7/-pairs. Possibilities of superconductivity due 
to the ?7-pairing mechanism were discussed in other related models, in particular, in one 
dimensional systems. See, for example, Refs. jH], and ^U] and references therein for more 
information. 

Another rigorous result is a theorem on the Hubbard model with attractive on-site in- 
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teractions on bipartite lattices due to S.-Q. Shen and Z.-M. Qiu [TT]. They proved that 
these models exhibit off-diagonal long-range order in the ground states. See also Ref. [12] 
for related results. 

At the present time, to draw an exact result about superconductivity in the Hubbard 
model with on-site repulsion or the t-J model seems to be a formidably difficult task, but we 
have a chance to obtain one in extended Hubbard models. Although these exactly solvable 
models have somewhat artificial aspects, they will give us some insight into mechanisms for 
the phenomenon. In this paper, motivated by this viewpoint, we propose a new exactly 
solvable model of strongly correlated electrons. It is based on a two dimensional d-p model, 
a tight binding model of electrons with on-site interactions on the decorated square lattice 
corresponding to the CuC>2 plane. The present new model has the infinitely large repulsive 
interactions on Cu-sites and furthermore contains correlated-hopping, pair-hopping and at- 
tractive charge-charge interactions of electrons. For even numbers of electrons we obtain 
the exact ground states of the model. A mechanism used to construct the ground states 
is similar to 77-pairing mechanism, but an electron-pair used here is different from 77-pair. 
Our electron-pair consists of singlet pairs of electrons on the same O-site, on the nearest- 
neighbor pair of the Cu-sites, and on the nearest-neighbor pair of the Cu- and O-sites. The 
constructed exact ground states are expressed as a linear combination of products of these 
singlet pairs and thus are regarded as RVB states. We also show that these ground states 
exhibit off-diagonal long-range order, which allow us to construct a ground state with explicit 
electron-number symmetry breaking. 

The recent experimental results of angle-resolved photoemission spectroscopy ^3] sug- 
gested that electron-phonon coupling is important in high-T c superconductors and thus a 
microscopic theory should include this effect. Although this fact does not immediately 
justify an addition of artificial interactions to a model Hamiltonian, we think that it is im- 
portant to know what kind of electron-electron interaction, whichever of Coulomb repulsion 
or electron-phonon coupling (or others) induces this, may stabilize a superconducting state 
in a theoretical point of view. We hope that our results will be useful in a future study. 

In the next section we define the model and state the main result as a theorem. In 
Section El we give some remarks, including the RVB representation and off-diagonal long- 
range order. In Section 0] we prove the theorem. 

2 Definition of the model and the main result 

We start by describing lattice A on which our model will be defined. Let L be an arbitrary 
positive integer and let D be a set of sites 

D = {x = (x (1) ,x (2) ) I x® G Z, 1 < x« < L for I = 1,2} (2.1) 
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with periodic boundary conditions. Let P be a set of sites located at the mid-points of 
nearest neighbor sites in D, 

p = { u = x + 5i/2 | x G D, I = 1,2}, (2.2) 

where 5i is the unit vector along the /-axis. Then we define A as A = D U P, which mimics 
the CuC>2 plane. For later use, we also define 

Px = {u | u G P, \u 

for x G D and 

D u = { x | x G D, \x 

for u G P. 

Let Q iCr and cj CT be the annihilation and the creation operators, respectively, for an electron 
with spin a at site i G A. These operators satisfy the usual anticommutation relations 



fa,*, c j:T } = {c[ a , c\ T } = (2.5) 

and 

{ c i,a, c), T } = bifior (2.6) 
for any i,j G A and a, r =|, |. The number operator n ij(T is defined as n it<T = c\ a c ija . We 



denote by N e the number of electrons in A and by <J>o a state with no electrons. 

We assume that repulsive interactions between electrons on D-sites are infinitely large 
and that each of these sites is at most singly occupied. So we consider the Hilbert space 
Hn s spanned by the linearly independent states of the form 




with arbitrary subsets A C D, B^,B± C P such that \A\ + \B^\ + \B^\ = N c . Here <r A is a 
shorthand for a spin configuration (cr x ) xeA . 

To define the Hamiltonian, we introduce the following new fermion operators: 

a u ,a = c U;a + a c x,a for u G P, (2.8) 

xeD u 

b x ,a = c Xj(7 + 13 °u,a for x G D, (2.9) 

ueP x 

d x ,a = E e- 2iQ 'X CT for (2-10) 



x| = 1/2} 
-u| = 1/2} 



(2.3) 
(2.4) 
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where a and (3 are real parameters with a ^ —j3, and Q = (0, tt). We also introduce operators 
n u a = a u a a u,a, n x a = o°x,o and n d x a = d x a d x , a - Then, by using these operators, we define 



H = V D [tY, E <* + * E E I Pd, 



(2.11) 



flint,! = -Vl E E <V»L " ^ 



E E E< 

ieDfr=T,J. xe-D 2=1,2 tr=T,i 



nx+5 Il( T + ^d, (2-12) 



and 



flint, 2 = v D -v 2 e E <-aV - w * E E E ^ (2-i3) 

\ ueP(T=T,| u£Px£D u a=i,i ] 

where t, s,Vi,Wi,V 2 and W 2 are real parameters, and —a denotes the spin opposite to a. 
The projection operator Vd which eliminates states with doubly occupied D-sites is defined 
by 

V D =X\V X (2.14) 



xeD 



with 



V x = (I- n x ^n Xt i). 
The Hamiltonian H is rewritten as 



(2.15) 



with 



H = V D [ E *<i c U c j,<r ] 



t if i = j e D; 

2pH + s if i = j e P; 

fit if|z-j| = §; 

(3 2 t if i j, i,j e P x . for some 

otherwise, 



(2.16) 



(2.17) 



(see Fig. ^) and it describes quantum mechanical motion of electrons feeling infinitely large 
repulsive interactions at D-sites. The Hamiltonians H mt) i and H mt) 2 correspond to correlated- 
hopping, pair-hopping, charge- charge, and spin-spin interactions. We can observe this by 
rewriting, for example, H 1X & 2 as 



fl 



int,2 



E ( - ft2 ^2 E E 4 -a^M,, 

\u€P x,y&D u ;xyty <r=1,l 



Figure 1: The lattice structure and the hopping matrix elements of H . The filled circles 
and the open circles correspond to P-sites and P-sites, respectively. The on-site potential 
of P-sites is t and that of P-sites is 2f3 2 t + s. The solid lines and the dashed lines represent 
the hopping matrix elements, f3t and /3 2 t, respectively. 

-aV 2 E (cx,-*Cu,-<r + c u ,-aCx ,-a)n u ,a - 2V 2 n u jn Ujl 
xeD u (j=T,i 

-±(a?V 2 + W 2 ) E n x n u + 2(a 2 V 2 -W 2 ) £ S< 3 >S< 3 >)W (2.18) 
xeD u xeD u I 

where rii = + and Si = (rijj — n^)/2. The Hamiltonian Pi n t,i can be rewritten as 
well, but that has a somewhat complicated form. We note that the interactions Pmt,i and 

i?int,2 conserve the electron number and the eigenvalue of J2ieA $i 3 \ but they do not possess 
the spin SU(2) symmetry. 

In this paper we consider the Hamiltonian 

H = H + Pf intjl + # int , 2 (2.19) 

with the parameters satisfying t — 4(1 + a 2 )Vi, W\ = a 2 V\ s = (1 + 2a 2 )V 2 , W 2 = a 2 V 2 and 
Vi, V 2 > 0. In this case, the direct spin-spin interaction terms vanish and the charge-charge 
interactions are attractive. 

In our construction of exact ground states, it is crucial to rewrite H as 

H=Y J H X +Y.H U (2.20) 

x£D u€P 
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with 

H x = V D (v, £ j V D (2.21) 

for x G D and 

H u = V D (v 2 cl a a u ^V D al_ a c u ,)j V D (2.22) 

for u6P. Here we used the operator identities Clj 

-V D c) a c i)(7 VD if i 7^ j; 

^d(1 - n iiT - n u )P D if z (2.23) 

[ PD(l-n ii(T )P D if i = j eP, 

and 



'PxCi,o'Px ^-i^cjPxi (2.24) 

V x c{ a V x = V x c{ a . (2.25) 

Operator identities (|2.24j) and ()2.25|) follow from straightforward calculations, and ()2.23|) is 
derived by using n x ^n x ^V x = and 

sp J (1 ^x,—u)c xcr 'P x {\ Tlx,— <j)Cx,<j Hi 2J) 2g") 

1 V x c^ a otherwise, 

-p t _ / 4, CT (1 - V<r) = c*,^ 1 - n x ,- a )V x if i = re; 

' " \ r! f ,P r otherwise, 1 j 

which also follow from simple calculations. We note that Hi with i e A are positive semidef- 
inite operators and thus a zero-energy state for all Hi (if it exists) is a ground state. 
Let us define 

C+ = £ e2iQ -"<T<l> ( 2 - 28 ) 

which corresponds to a singlet two-electron state. The main result in this paper is the 
following theorem: 

Theorem 2.1 Fix the electron number N e and suppose that N c < 2\P\. When N e is even, 
the ground state <&g,n c °f H is unique and is given by 

®g,n c = V D (C f ) 2 $0 = (VdC) 2 $0, (2.29) 

which satisfies H^q ^ c = 0. For odd N e , the ground state energy is positive. 

It is noted that the ground state energy for odd N e may converge to zero as L — > oo. 
Whether this is the case or not should be clarified in a future study. 
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3 Some Remarks 



It is easy to verify the relation 

C f = E E ^ ("V^, + afl +Si/2 + afl_ Si/2 ) + E e 2iQ - u < T ci a , (3.1) 
xeD 1=1,2 ueP 

where fjj = cj jcj j + c^c]^ is the creation operator for the valence-bond (singlet pair) on 
sites i and j. Furthermore the projection operator Vd eliminates the terms of the form 

fl,iflj---*o = -<iAifli-'-*<>> ( 3 - 2 ) 

when we express &g,n c by using f-j and cj^cj^. Thus the ground state $g,jv 6 has the form 
of the RVB state which is a linear combination of products of on-site singlet pairs cj^c^ 
in P and the valence-bonds fl^ x+5l and fl x±S[ / 2 - ft is noted that the weight for the singlet 
pairs fl x+ s, on the nearest-neighbor D-sites becomes relatively large for 1/a <C 1. 
The term —V\n x _ a n xa in interaction Hamiltonian ifi n t,i is rewritten as 

" VX,-X,a = -4(1 + « 2 )(1 + 4/3 2 )^ + 4(1 + a^V^bla + V 1 d x , JU'r.A,, a- (3.3) 

From this expression, we find that — V\n d x _ a n x a is bounded below by —4(1 + a 2 )(l + 4(3 2 )Vi 

and furthermore that it is attained by states of the form b x c d x _ c • • • <E»o- Therefore this 
interaction term is describing an attractive and magnetic interaction between localized single- 
electron states corresponding to fermion operators b x>a and d XjlJ . The other terms can be 
rewritten into similar forms and interpreted as attractive and magnetic interactions as well. 
Our model says that attractive magnetic interactions between localized states together with 
the strong on-site repulsion can stabilize the condensed RVB states consisting of local singlet 
pairs. 

The explicit expression of the ground states makes it possible to evaluate order parame- 
ters. We first note that, in the case of a — 0, $g,tv c i s reduced to the ^-pairing state on the 
P-sites, and therefore the ground states have off-diagonal long-range order. 

For a^Owe can prove the existence of off-diagonal long-range order as follows. Let us 
introduce order parameters 



and 
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Noting the anticommutation relation 

{c u ,a, a u >,<j} = Ku' (3.6) 
for u, u' G P, one finds the commutation relation 

[A 1)A , VdCVd] = V D El 1 - <T C «>T - <|^l)J ^ (3-7) 

and 

V D £ (1 - < tCu , t - 4 )iCu4 )^$ GliVe = (|P| - iV c )$ G ,iv c . (3.8) 
These two relations lead to 

Al,A$G,Ar e = - f - l)) $G,^-2. (3.9) 

Therefore we obtain 

(^^(y-ByIy- 1 ))' < 3 - 10 ' 

where the expectation value (• ■ -)a,jv i s defined by 

/ v ($G,7V e , • • • $G,V C ) /o ii\ 

(• • ■)A,N e = —7^ = (3.11) 

(.*G,JV e) ®G,N C ) 

We write [ia,n c for the right-hand-side of ()3.10|) . Then, by using the Schwarz inequality 

r < (a^ a a 2 

AM.iVe A)A,7V e > (3-12) 
and inequalities (A{ A Ai >A ) Aj jv e < 1 an d (A^ A A2 )A ) Aj jv e < (1 + 2a) 2 , we find 

"~ VX < (A{ A A hA ) A , Ne < 1 (3.13) 



(l + 2a 2 ) 2 
and 

A, Ne < (A5, A A 2)A ) AiiVe < (1 + 2a 2 ) 2 . (3.14) 

Let {$Gi}i =0 be a sequence of ground states &g,n c on the lattice with side-length L such 
that the electron filling N e / (2|A|) = iV c / (6L 2 ) converges to v as L — ► oo, and let {(• ■ -)^}iLo 
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be a sequence of expectation values with respect to $q £ . Then from inequalities ()3.13j) and 
(|3~T3J) we have 

WW £ (TTW (3 - 15) 

and 

liminf(A^ A A 2iA )(r ) > / i2 (3.16) 

where /i„ = ^ (l — if) > w ^ich imply the existence of off-diagonal long-range order for < 
v < 2/3. 

Now we can construct a ground state with explicit electron-number symmetry break- 
ing [T3]. Let us define 

A t 

&G,N. = ®G,N C + j 2A $G,7V C , (3.17) 

V ( A 2,AA 2)A ) Ai Ar c 

which is also a zero-energy state. We note that the limit infimum of (A 2 ,aA 2 A }^ ^ is bounded 
below by nl, since 

(A liA A 2iA } A)JVe = (A| !A Ai )A ) A) jv e + i^pd^l " ^e), (3.18) 

which follows from f)3.7|) . Since two states with different electron number are orthogonal, we 
find that 

( $G,iV c! A2,A $ G ,aQ _ 1 
($G,iV e , $G,JV ( , 

Therefore, the limit infimum of the sequence of (A 2jA ) A n obtained by using {^glIl^o w hh 
f!3.17|) is bounded below by fi u /2, which implies electron-number symmetry breaking. 

We finally remark about extensions of the present model. It is possible to construct 
similar models in one, three and more dimensions by the same method. The model in three 
or more dimensions may exhibit a finite temperature phase transition. It is also possible 

to construct models whose ground states are written as 2 $ . The details will appear 
elsewhere. 

4 Proof 

Proof of Theorem \2.1\ We first prove that 3>g,n c m ()2.29j) is a zero-energy state for all Hi 
and thus is a ground state. Using the anticommutation relations (|3.6J) and 

{cx,a,al } = ax[u G P x \ (4.1) 



(a 2 ,a>u = v t, T 7 = ~ v^-a^aK^- (3.19) 
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for x G D, where x[" even t"] takes 1 if "event" is true and takes otherwise, one finds the 
following two commutation relations: 

[c u , a ,^}=^ u ai,-, (4.2) 

for u G P and 

[ CavT) C f ] = aa E e^al (4.3) 

for x E D. (In the right-hand-sides of (|4.2jl and (|4.Hjl . the coefficients a =f and J, are 
regarded as +1 and — 1, respectively. We will use this convention in the following.) By using 
operator identities (|2.25|) . (|2.26jl and commutation relation (|4.2j) . we find 

{V D a ] u _ a c u ^)V D C = V D ai_ a c u ^C ] 

= V D al_ a (^c u , a + ae 2iQ - u al_ a ) 

= V D C,\V D a{^c u ^. (4.4) 

To get the third line, we used (a^ ff ) 2 = 0. This implies that (PL>al,- o -c U;CT )(P£)C^)^ < ^0 — 0, 
and thus H u $ G Nc = for all u G P. Similarly, by using (|2.25|) . (|2.26|) and commutation 
relations (JOJ, (Oil we have 



V D di_ a E c u , a V D C = V D dl_ a E c u,. C 1 " 
ueP x J V " eP * / 

u£P x ueP x 

= V D C \V D dl_ a E c u A (4.5) 

V udP x j 

and 

= V D {1 - n x _ a )d ] x _ a ( ] Cx,a 

+aaV D (l - n x ^)dl_ a E e 2iQ M,-* 

ueP x 

= V D (1 - n x ,. a )C (V D d{_ a c x , a ) . (4.6) 

Here we also used {c x - a , d], „} = and {d x _ a ) 2 = 0. It follows from ()4.5|) and ()4.6|) that 

(P^c?J. j _ .6 Xi(J )('PdC^)^ ( ^o = 0, i.e., if x .$G,7v c = for all x G D. Therefore we conclude that 
H$ G , Nc = 0. 

The proof for the other statements in the theorem relies on the following lemma, which 
will be proved later. 
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Lemma 4.1 Any zero-energy state $ for H in the Hilbert space 07v o = o Hn c is written as 

e *wpd ( n <t) f n $ o, (4.7) 

BCP VmSB / \ueB ) 

where <f){B) are real coefficients. Furthermore coefficients satisfy \<fi(B)\ = \<f)(B')\ for any 
B,B' C P such that \B\ = \B'\. 

It immediately follows from Lemma 14. II that the ground state energy for odd N e is always 
positive. The remaining task is to prove the uniqueness. Suppose that there are two zero- 
energy states for fixed even N e . Then, an arbitrary linear combination of these states is also a 
zero-energy state, which should satisfy Lemma f4. 11 However, we can make a suitable linear 
combination so that a coefficient <P{Bq) for a subset Bq will be vanishing, and this leads 
to the conclusion that all the other coefficients are also vanishing. This is contradicting 
with the assumption, and therefore the ground state is unique. This completes the proof of 
Theorem 12.11 ■ 

Before proceeding to the proof of Lemma 14.11 we prove the following lemma. 
Lemma 4.2 For A C D, spin configuration cr A and B^, Bi C P, define 




Then, the states Vd&(A, <t a\ B^, B^) are linearly independent, and the collection of these 
states with \A\ + \B^\ + \Bi\ = N c spans the Hilbert space Hn c . 

Proof of Lemma \4.^\ Consider a linear combination 

$ =EE E ^A^B^Vd^A^B^B^ (4.9) 

AcD (T a B ] ,B 1 CP 

where J2cr A means the sum taken over all spin configurations. We suppose that $ = 0. From 
the anticommutation relation (jH.fi)) it follows that the inner product 

(*(A', (t' a ,- B' p B{), V d $(A, ct a - B h B{j) (4.10) 

is zero if both Bi C B^ and B[ C B^ do not hold. (See (|2.7j) for the definition of 
ty(A, cr A ; B^, Bi).) We furthermore have that 

{9 (A', er' A ,; £? T , B{), V D ®{A,* A] B h B l ))=U (4.11) 

for A' ^ A and that 

(*(A, <r' A , B h V D $(A, a A ; B h BJ) = X [< = °x for all x 6 A]. (4.12) 
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Thus we obtain from fit {A, er A ; P,P), $) = that </>(A, cr A ; P, P) = for any A C D and 
spin configuration cr A . Then, examining 

(*(A,<t a ;P\{u},P), $) = (4.13) 

for « 6 P, we find that 4>(A, cr A ; P\{u}, P) = 0. Similarly, examining 

ct a - P\B h P\Bi), $) = (4.14) 

for 5|, S| C P with |JB||, = 0, 1, 2, . . . |P| repeatedly, we conclude that all the coefficients 
<r A ; -Bf, Pj) are vanishing. This proves the first claim, and the second claim is now 
trivial. ■ 

Proof of Lemma \4-l\ Suppose that $ is an arbitrary zero-energy state for H in 0jy^l o H/v c - 
By using the basis states Vd®(A, er A ; B^, B±) with \A\ + \B^\ + |£>J < 2\P\, we represent $ 
as 

$ =EE E 0(A^;P T >^)^(A^;P T >^)- (4.15) 

The sum is restricted to |v4| + |P|| + |PJ < 2\P\, but, here and in the following, we do not 
written down this restriction for notational simplicity. 

Since the local Hamiltonian Hi is the sum of two positive semidefinite operators, zero- 
energy state $ for H must be annihilated by these operators for any i G A. 

We first examine the case i = u G P. In this case the condition H u § = is equivalent to 

V D al_ a c u , a V D $ = (4.16) 
for a =T, !• For a =f the left hand side of (|4.16|) becomes 

X0(A, <t a ; P t , P x )P d $(A <x a ; fl T \{u}, #1 U {«}), (4.17) 

where sgn[w; BJ is a sign factor arising from exchanges of fermion operators. Since all 
the terms in (|4.17|) are linearly independent, we have the condition 4>(A, a A \ B^, B\) = for 
P|, Bi such that u G B<\ and u ^ Pj. A similar calculation for cr = j yields cr A ; Bf , B\) = 
for B^,Bi such that u ^ B^ and w G -B;. Since these conditions must be satisfied for all 
u G P, we obtain <p(A, <r A \ B^, B±) = for B^ ^ B^. 

So far we have shown that a zero-energy state $ can be expanded as 

$ = E E E <KA p)Pd$(A o-a; a, b), (4.18) 

AcD CT A _BcP 
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where 4>(A, cr a\ B) = 4>(A,<ta', B, B). We further derive conditions on <f>(A, cr a\ B) from 
Hi<& = with % = x G D, which is equivalent to 

V D dl_ a b x ,aV D ^ = (4.19) 

for cr =|, |. The left-hand-side of ()4.19j) can be expanded by using ()4.8j) . so that we decom- 
pose this as 

E E E •••+ E E E ••• (4-20) 

A'cD;x<£A' <T' A , B^B^cP A'cD;x&A' (T' A , B^B^cP 

and write $ x for the first sum. Since the states in the first sum and those in the second sum 
are linearly independent, 3> x = must be satisfied. 

To obtain Q x for a =f, we operate Vd(^ x ^x,{Pd on the basis state Vd^{A,cta', B, B). 
Then, we have 

(l-n xd )V D dl l c x ^(A,a A ;B,B)+pV D dl l I ]T c u A $(A,* A ;B,B). (4.21) 

\UGP X J 

The second term in (j4.21j) becomes 

(3 Y, sgn[u,u';B}x[ueB]x[u' <£ B]e 2i Q u 'V D $(A, <r A ; B\{u}, B U {u'}), (4.22) 

u,u'eP x 

where sgn[w, u';B] is a sign factor arising from exchanges of fermion operators. Thus, the 
second term contributes to <& x only when x A. The first term in ()4.21|) becomes 

(X[x E A] X [a x =T] + x[x i A])V D dl lCx ^{A, a A ; B, B) 

-(X[x e A)x[a x =T] + x[x i A\)V D d\ A c\ A c xX c x ^{A, <r A ; B, B). (4.23) 

Since site x is always occupied by an electron in the second term in ()4.23|) . this term never 
contributes to $ x . The first term in ()4.23j) furthermore becomes 

X[x G A]x[(Tx =T]sgn[ar;A] ^ sgn[w'; B]x[u' £ B] 

u'eP x 

xe 2iQ u 'V D $(A\{x}, (TA\ {x y, B,BU {«'}) 




+X[x £ A]a s S n h u'\ B\ X [u G B] X [u' $ B] 

u,u'eP x 

xe 2iQ - u 'V D $(A, a A ; B\{u}, B U {«'}), (4.24) 
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where sgn[x; A] and sgnfw; B] are again sign factors. Since site x is always occupied by an 
electron in the second term in the above expression, this term does not contribute to $ x . 
Therefore we finally obtain 

$ * = E E E sgn[x;A] sgn[u';B}x[u' <£ B] 

AcD;x&A(T A ;a- x =] BcP u'£P x 

xe 2[Q - u '<P(A, <r A ; B)V D <5>(A\{x}, ct A \ {x} ; B,BU {«'}) 
+ E EE E s § n K u- B] X [u G B] X [u' i B] 

AcD;x(£A (J a BCP u,u'£P x 

xe 2iQ -"'0(A (T A ; B)V D $(A, <r A ; B\{u}, B U {«'}). (4.25) 

It is noted that the terms in the first sum are linearly independent of those in the second 
sum. 

Choose a configuration (A, a A \ B) satisfying that A contains x, a x in a A is j and there 
exists u' G P x such that v! $ B. Then, by checking the coefficient of the basis state 
Vd&(A\{x}, <T J 4\{ a! }) B, BU{u'}) in (|4.25|l . we obtain 4>(A, <r A) B) = for such configuration. 
Since this and a similar result for a —I must hold for any x G D, we obtain that 

<P(A,<r A ;B)=0 (4.26) 

if (U xeA P x ) is not a subset of B. 

Next, choose a configuration satisfying that A does not contain x and there exist u, v! G P x 
such that u G B and u' £ B. For this kind of configuration, by checking the coefficient of 
the basis state P D $(A, cr A ; B\{u}, B U {«'}) in ()4.25j) . we have that 

(a + /3)(sgn[u, B]e 2iQV 0(A <r A ; B) 

+sgnK, u; 5^^]e 2iQ ^(A tr A] B u ^ u ,)) = 0, (4.27) 

where B u ^ u i = (B\{u}) U {u'}. Since sgn[-], e 2l< ^"" and e 2l< ^""' give only sign factors and a + /3 
is non-zero by definition, we obtain the (necessary) condition that, for any A not containing 
x and any cr A , 

\<j>(A, a A - B) | = \(j)(A, ct a - B u ^ u ,) | (4.28) 

if there exist u, u' G P x such that u G B and v! B. 

By using conditions ()4.26|) and ([4.28)1 . we shall prove that 4>(A, a A \ B) = for any A ^ 0. 
Suppose that <f>(A, er A ; B) ^ for some non-empty set A C D, spin configuration cr A and 
BcP. Here we say that w and w' in B are connected when D u fl 7^ and decompose 
B into connected components as B = B 1 U • • • -B m . We find from condition (|4.2fij) that A is 
contained in U ug £.D M , and without loss of generality, we can assume that U ue B 1 P J u contains 
at least one site in A. For connected subset B 1 we can find at least one pair of sites uq, U\ G P 
and a site x\ G U u( zb±D u such that uq B l ,u\ G B l and Uq,ui G P Xl . When m > 1 this 
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is trivial, and when m = 1 this follows from N e = \A\ +2\B\ < 2\P\. It is noted that x± is 
not in A due to condition ()4.26|) . Then, since we can always find subsets {u 2 , ■ ■ ■ , u n } C B 1 
and {x 2 , ...,x n }c U ueB i D u satisfying that u\ (fc B^ l _ yUQj Ui+i G B l Ui ^ U0 , u h ui +1 G P Xl+1 and 
xi A for 1 < I < n, and x n G A, the repeated use of (|4.28jl gives 

\(f>(A,tr A ;B)\ = tr A - B Un _^ Uo )\. (4.29) 

But, since u n -\ G P Xn is not in _B Un _ 1 ^ U0 , the coefficient (f>(A, cr A \ B Un _ 1 ^ U0 ) must be vanish- 
ing because of (|4.2fij) . which leads to a contradiction. Therefore the claim is proved. 
From the above results we conclude that $ is expanded as 

$ = J2 4>(B)V D $((b, <r ; B, B), (4.30) 

BCP 

where <fi(B) = 0(0, erg; B). Then, using ()4.28|) again, we find that |0(-B)| = \<p(B')\ whenever 
\B\ = \B'\, which completes the proof of Lemma f4. II ■ 
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